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Abstract 

We consider a quasilinear parabolic differential equation associated with the renormaliza- 
tion group transformation of the two-dimensional hierarchical Coulomb system in the limit as 
the size of the block L J. 1. We show that the initial value problem is well defined in a suitable 
function space and the solution converges, as t —* oo, to one of the countably infinite equi- 
librium solutions. The j-th nontrivial equilibrium solution bifurcates from the trivial one at 
(3j = 8ir/j 2 , j = 1, 2, . . . . These solutions are fully described and we provide a complete analy- 
sis of their local and global stability for all values of inverse temperature > 0. Gallavotti and 
Nicolo's conjecture on infinite sequence of "phases transitions" is also addressed. Our results 
rule out an intermediate phase between the plasma and the Kosterlitz-Thouless phases, at 
least in the hierarchical model we consider. 



1 Introduction 

We consider, for each (3 > 0, the partial differential equation 



u t - -^-(u xx - u 2 x ) - 2u = (f.f) 

47T 

on R + x (— 7r,7r) with periodic boundary condition, u(t,—7r) = u(t,ir) and u x (t, — n) = u x (t,n), 
in the space of even functions, satisfying an additional condition u(t, 0) = 0[|. We show that the 
initial value problem is well defined in an appropriate function space B and the solution exists 
and is unique for all t > 0. Furthermore, as t —>■ oo, the solution converges in B to one of the 
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1 This is assured by a Lagrange multiplier (see Remark 3.1 
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(tt) and 0'(-tt) = 0'(tt). For /? > 8tt, ■ 
. For ft < 8n such that 8n / (k + l) 2 



6 = is the (globally) asymptotically stable 
< ft < 8n/k 2 holds for some k e N A 



unstable and there exist 2k non-trivial equilibria solutions 
the only asymptotically stable ones. 



<+, <?>o is 
± 



of ( |1.2| ) among which 0^ are 



^ have a (j — 1) -dimensional unstable 
As a consequence, there exists a dense 



The aim of the present work is to show that, for j > 1, 
manifold Aij C £> so 4>f are more stable than 0^ if j < j'. 
open set of initial conditions in B such that (fif (0]~ is not physically admissible) is the non-trivial 
stable solution for all ft < 8n. 

Our description of equation (|1.1|) is motivated by two distinct goals. Firstly, it provides a new 
example of nonlinear parabolic differential equation by which a geometric theory can be carried 
out (see e.g. Henry [[H]]). According to this theory, the above scenario can be stated as follows: 
there exist a sufficient large ball B C B about the origin such that, if u(t,B ) denotes the set of 
points reached at time t starting from any initial function in Bq, then the invariant set f] t>0 u(t, Bq) 
coincides with the fc-dimensional unstable manifold JCk = Uo<j<fc ^-i = provided 8n / (k+1) 2 < 
ft < 8n/k 2 . 

Secondly, the solution of the initial value problem ( |1 . 1| ) describes the renormalization group 
(RG) flow of the effective potential in the two-dimensional hierarchical Coulomb system and the 
stationary solutions {0^}, the fixed points of RG, contain informations on its critical phenomena. 

The analysis of equation ( |1.1[ ) presented here can hopefully bring some light to a question raised 
by Gallavotti and Nicolo [|GN|| on the "screening phase transitions" in two-dimensional Coulomb 
systems. The existence of infinitely many thresholds of "instabilities" found in the Mayer series at 
inverse temperature ft n = 8n(l — l/(2n)), n G N+, indicates, according to the authors, a sequence 
of "intermediate" phase transitions from the plasma phase (ft < ft\ = An) to the multipole phase 
(ft ^ ftoo — 8 71 ")- They conjectured that some partial screening takes place when the inverse 
temperature decreases from 8n to 47r, which prevents the formation of neutral multipole of order 
larger than In where n is the integer part of 1/(2 — ft /Air) (dipoles are the last to be prevented at 
An). 

The Kosterlitz-Thouless phase (multipole phase) was established by Frohlich-Spencer ||FS|| and 
extended up to 8n by one of the present authors and A. Klein |[MK|| . Debye screening (plasma 
phase) was only proved for sufficiently small ft << An |[BF|| . Study of the region [An, 8n] began 
with the work by Benfatto, Gallavotti and Nicolo BGN | on the ultraviolet collapses of neutral 
clusters in the Yukawa gas which served as a base for the results in [GIN]. It seems improbable, 
on the light of the present knowledge, that a conclusive answer to the Gallavotti-Nicolo conjecture 
will come up soon. It may be noted, however, that the scenario of an intermediate phase, which 
has challenged the conventional picture due to Jose et al ||JKKN 
et al [ FLL ] based on Debye-Huckel-Bjerrum theory and by Dir 
reinterpreted the ultraviolet collapses in the Yukawa gas. 

The Kosterlitz-Thouless phase is manifested in the hierarchical model as a bifurcation from the 
flMFH 



has been contested by Fisher 
who have 



trivial solution 



Our results rule out the existence of further phase transitions since no other 



bifurcation arises from the stable solution (see Theorem |5.1| on the stability of <pf). 
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Even though the existence of the invariant unstable manifold ICk may provide a suitable expla- 
nation to the appearance of Gallavotti-Nicolo's thresholds, the nature (and location) of the insta- 
bilities in the hierarchical Coulomb gas differs substantially from the one we have just described, 
because neutral multipoles cannot be formed in the hierarchical model. We believe, however, our 
investigation may be helpful for the plasma phase. Numerical analysis shows the stable solution <p± 
looks like the Debye-Hiickel potential <Pdh = (2vr//3)a; 2 in (— 71,71) right after the transition takes 



place (see Remark fj~6| ) . 

As in |F), the renormalization group (GR) flow ( |1 . 1|) may be derived from the block-spin RG 
transformation of a two-dimensional hierarchical Coulomb system in the limit as the block size 
L I 1. This procedure, called local potential approximation, has been discussed by Felder ]F] in the 
context of Dyson's hierarchical model, whose partial differential equation, 

1 d-2 1 2 , , 

u t --u xx ^ — xu x - du + -u x = , (1.3) 



coincides with (|l.lj) when his dimensional parameter d = 2 if [3 is equal to 2tt (without boundary 



conditions). Felder showed that ( |1.3|) has global stationary solutions u\ n on R for 2 < d < d n with 
u 2n( x ) ~~ * as d T d n and calculated their profile. Here, d n = 2 + 2/{n — 1), n — 2, 3, . . . , is the 
sequence of thresholds where nontrivial fixed points are expected to appear as a bifurcation from 
the trivial solution. We mean by global solution one which doesn't blow up at finite x. 

The present paper begins with a derivation of equation (1.1) in Section [| The existence, 



uniqueness and continuous dependence on the initial value are presented in Section R] and the 



precise statements are given in Theorems |3.2| and |3]4]. We describe all global solutions of ( |1 . 2| ) 
completely in Section |j. Due to smoothness and the periodic condition, blow-up of an admissible 
stationary solution is impossible. We show that the non-trivial stationary solution for /5 < 87r is 
unique modulo solutions with period 2n/j, j = 2, 3, . . . , which are responsible for the existence of 
the unstable manifold (see Theorem (4.1|) . Finally, we analyze in detail the local and global stability 
of equilibrium solutions of (|1 . 1| ) in Section |5]. The main results are stated in Theorems |5J] and 



2 The Flow Equation 

This section is devoted to the derivation of ( |1.1| ) from the RG transformation of two-dimensional 
hierarchical Coulomb system. We begin with a brief review of this model. 

A Coulomb system is an ensemble of two species (for simplicity) of charged particles, interacting 
via a two-body Coulomb potential V. In the grand canonical ensemble the total number of particles 
fluctuates around a mean value determined by the particle activity z. It will become clear that the 
charge ensemble, rather than the particle ensemble, is more appropriate for RG transformation. 

A configuration q of this system is a function q : A C Z 2 — > Z which associates to each site x 
of the lattice A the total charge q(x) at this position. 

To each configuration we introduce two functionals: the total energy E : Z A — ► M + , 

E(q) = \ J2 <l(x)V(x,y)q(y) (2.1) 

x,yeA 
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(self-energy is included) and an "a priori" weight F : Z / 



F(q) = l[X( q (x)) 



(2.2) 



defined for positive real valued functions A. 
The equilibrium Gibbs measure /ia : Z A 



is thus given by 



1 



HA(q) := —F(q) e 

—A 



-/3B(?) 



where /3 is the inverse temperature and 



-A 



-PE(q) 



[2.3) 



(2.4) 



is the grand partition function. 

It has been shown (see e.g. ||FS|| ) that the standard Coulomb system in the grand canonical 
ensemble with particle activity z has charge activity given by X(q) = I q (2z), where I q is the g-th 
modified Bessel function. If X(q) = 6 qt o + z (5 q ^i + S q -x), 5a is the grand canonical ensemble of 
charged particles with hard core. 

Let us introduce our hierarchical model as proposed in ref. ||MP| . The potential V in ( |2 . 1| ) is 
replaced by a function 

V h (x,y) = -—lnd h (x,y) , 

given by the asymptotic behavior of the two-dimensional Coulomb potential with the Euclidean 
distance \x — y\ replaced by hierarchical distance 



d h (x,y) :=L N ^ 

defined for an integer L > 1, where 

N(x,y) := ini\N G N + 



x 



L N 



L N 



(2-5) 



(2-6) 



and [z] G Z 2 has components the integer part of the components of z 6 R 2 . Notice that dh is not 
invariant by translations. 

Now, given an integer number N > 1 , let A = A^r = [— L , L N — L^ -1 ] 2 fl Z 2 and define, for 
each configuration q G Z A , the block configuration q 1 : Ajv-i — ► Z, 



= i( Lx + y)- 



(2.7) 



0<Vi<L 

i=l,2 



The renormalization group transformation TZ acting on the space of Gibbs measures 



3D, 



f2.8) 



N ; 

fixed 
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involves an integration over the fluctuations about q 1 following by a rescaling back to the original 
lattice. 



As it has been shown in ||MP|| , the RG transformation 1Z preserves the form of the Gibbs measure 
in the grand canonical ensemble of charges. The measure fi\ is thus given by (|2.3| ) with the "a 
priori weight" F replaced by 

F 1 (q 1 )= II A V(*)) ( 2 - 9 ) 

xeAjv-i 

where 



A 1 (p) = L-* 2/(47r) j\ ★ A ★ • • • ★ A ) (p) (2.10) 

L 2 —times 

with (A -k g){p) = A(p - O) q{q)- Note that H Ajv (A) = S Ajv _ 1 (A 1 ). 



Remark 2.1 A peculiar feature of hierarchical models is the reduction of the measure space where 
TZ acts to local functions. The RG transformation j\2.%) induces a transformation A 1 = rA given by 
(\2. 1 (\) on the space of infinite sequences. Note that the space £i(Z) of summable sequences is closed 



by the r transformation: (A * A) G ^i(Z) if A G £i(Z) i/ie Hausdorff- Young inequality. 

In order to take L J, 1 limit of the RG transformation r it is convenient to write the system in 
the sine-Gordon representation. Fourier transforming (|2.10|) , 

%) = EA(?)e^, 

and using the convolution theorem, yields 

A%) = r\((p) = ±- f *(<p ~ r) \ L \r) dr (2.11) 

where 



1 _ ^ c -7r(y+27r») 2 /(/31nL) ^.12) 



(/JlnL) 1 ^ 
by the Poisson formula. 

Plugging Q2.12D into ( 2.11|) and changing the variable ( = r + 27m, equation ( j2.11j) can be 
written as 



r\(cp) — ( 1/ * A (» (2.13) 
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where u* means convolution by a Gaussian measure with mean zero and variance /31nL/(27r): 



-1/2 



/(C) 



(2.14) 



= e (/31nL/47r)(d 2 /V) 

where in the second form of the Gaussian convolution we have used Wick's theorem. 



Note that (|2.13| ) is precisely the RG transformation derived by Gallavotti who has started 
directly from the sine-Gordon representation. 

In order to let the block size L to 1, we introduce a variable t := n InL which keeps track of 
the number of times the RG transformation ( |2.8| ) has to be iterated in order to bring two sites at 
hierarchical distance L n to 0(1) distance. We shall take the limit L j 1 together with n — > oo 
maintaining t fixed. 

Define 



u(t, x) — — ln X n (x) 



(2.15) 



where \ n = r n \ denotes the n -th iteration of the transformation Q2.13 ). If one writes t' 
(n + 1) InL then, by taking the logarithm and using (|2.15|) , equation ( p,13|) reads 

W n , f ( Pt d 2 

u(t , x) 



u 



u 



which, combined with 



1 "{^(^^) exp( - e2,/ "" (f ' l)) 

{t, x) - ln |l + l - (u 2 x (t, x) - u xx (t, x)) - 2u(t, x)j + O (J^J | (2.16) 
(t, x )+ t -(j^ { u ^(t, x) - u 2 x (t, x)) + 2u(t, xfj + O (J-^J 

u(t', x) — u(t, x) 



u t (t,x) = lim 

v ' ; fit t'-t 



lim — (u(t r , x) — u(t, x)) 

n— >oo t 



(2.17) 



yields equation (|1 . 1|) . 



3 Existence, Uniqueness and Continuous Dependence 

In this section the existence, uniqueness and continuous dependence on the initial value of equation 
(|1 . 1|) will be established by Picard's theorem for Banach spaces. 

To avoid the appearance of zero modes upon linearization, we differentiate ( |1 . 1| ) with respect 
to x and consider the equation for v = u x , 

v t - j- (v xx - 2vv x ) - 2v = , (3.1) 
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with v (t, — tt) = v (t, 7r) and v x (t, —it) = v x (t, tt), in the subspace of odd functions and initial value 
v(0, •) = v . Note that the operator defined by the 1. h. s. of ( ^.1| ) preserves this subspace. 

Before we proceed, we have the following 

Remark 3.1 The "a priori weight" \{t,q) := A™(g) at scale t = nlnL, is a positive symmetric, 
X(t, q) = \(t, —q), sequence of real numbers and has to be normalized at all scales. In \MF\j equation 
t \2.1Q) was redefined so that A™(0) = 1 holds for alln. Here, the appropriated normalization is given 
by 

$>(t, g ) = i, 



since, in view of equation $2.13i ), this leads to the condition u (t, 0) = 0, which is already imposed 
for all t if 



(3.2) 



u(t,x) = / v(t,y) dy 
Jo 

with v(s,x) an odd solution of From we have 

u t = v t (t,y)dy 



(3.3) 



[a (u xx - u 2 x ) + 2u] x dy 

= a (u xx - u 2 x ) +2u- au xx (t, 0) 

where u x (t,0) = v(t,0) = by parity. Note that u(t,x) = — lnA n (x) + lnA n (0) also satisfies $31% ) 
by equations ($.ldj ) and ( $.17[ ). Moreover, note that there is a one-to-one correspondence between 
the solution of ( \1.1\ ) and the solution of (3J5), with the same initial value uq, given by 



and 



u(t, x) = u(t, x) — u(t, 0) 



u(t, x) = u(t, x) + a / e 2< -* s ^u xx (s, 0) ds 



(3.4) 



(3.5) 



where au xx (t,0) is the required Lagrange multiplier introduced in (3~5) to assure that u(t,0) = 
(see comments after equation (1.1') in ref. fiFtf). This correspondence will be useful in Section 

Because the standard initial condition uq(x) — z (1 — cosx) satisfies u' (0) = u' (tt) = 0, equa- 
tion \3. 1]) may equivalently be considered on (0, tt) with Dirichlet boundary conditions v (t, 0) = 
v(t,ir) = 0. 

Another reason for considering ( [3.1|) instead of ( |1.1|) is the fact that the nonlinearity 2vv x is 
more suitable than u 2 x for the analysis of equilibrium solutions and corresponding stabilities given 
in the next sections. 
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The boundary and initial value problem (|3.1|) may be written as an ordinary differential equation 

dz 

f t +Az = F(z) (3.6) 

in a conveniently defined Banach space B where 

Az = -az" - 2z and F(z) = -2az'z, (3.7) 

with a = (3/(4tt) and initial value z(0) = z$. 

The linear operator A is defined on the space C 2 of smooth odd and periodic real-valued 

f(x) g(x) dx. Because of (/, Ag) = (Af, g) , 

A may be extended to a self-adjoint operator in L 2 p (— vr, 7r). The domain D(A) of A is 

D(A) = {fe L\ v (-vr, 7r) : Af E L 2 p (-vr, tt) } 

and the spectrum of A, 

cr(A) = {X n = an 2 -2,nG N + } , (3.8) 

consists of simple eigenvalues with corresponding eigenfunctions (f) n (x) = (1/tt) 1 ^ 2 sin nx . 

Let Ai denote a positive definite linear operator given by A if a > 2 and A + a/ for some 
a > 2 — a, otherwise. The following properties also hold for A given by the closure in L q o p (— n, n), 
1 < q < oo, of the operator (—ad 2 /dx 2 — 2)| c2 



-o,p 



1. The operator A generates an analytic semi-group T(t) = e * A given by the formula 

T(t) = — [ —^—e xt d\ 

where T is a contour in the resolvent set of A with argA — > ±6, 7r/2<^<7r, as |A| — > oo. 
From this, we have 

\\e- tA \\ < C e~ ct and \\Ae- tA \\ < - e~ ct (3.9) 



t 



for t > 0, c < infA a (A) and C < oo. 
2. Given 7 > 0, let the fractional power of A\ be given by 

1 r°° 

A-~< = ^- / V- l e~ Mt dt 



r(7) 7o 

and define AJ = (A^ 1 ) \ A 1 7 is a bounded operator (compact if 7 > 0) with v4 x (d/ dx) 

and (d/dx) A l bounded in the L„ p (— 7r, 7r) norm. In addition, for 7 > 0, A 7 is closely 
defined with the inclusion D{A\) C -D(AJ) if 7 > r. 



2 From here on, the subindexes in Cq P , £o P , £c P , fle P and e<c, indicate spaces of odd and periodic (o,p) or 
even and periodic (e,p) functions. 
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It thus follows from 1. and 2. (see e.g. Q) 

\\Ale- tAl \\<^e- ct (3.10) 

holds for < 7 < 1, t > 0. Here C 7 is bounded in any compact interval of (0, 1) and also bounded 
as 7 \ 0. Note that, if the operator norm is induced by the L 2 -norm, equation ( 3.10 ) hold with 



C 7 = sup |(tA n ) 7 e~ tXn \ < sup |r 7 e- r | < (-Y , (3.11) 

nGN + r>tc \&/ 

uniformly in 7, t > 0. 

Following Picard's method, let us replace F in Q3.6|) by a locally Holder continuous function 
/: 0.V B: 

\\f(r)-f(s)\\<C\r-s\ e 



for < r < s < T and 9 > 0. In this case, a solution to ( |3.6|) is given by the variation of constants 
formula 

z (t) = e~ tA z + [ e~ {t ~ s)A f(s)ds. (3. 12) 

Jo 

Note that z : [0,T) — > i3 is continuously differentiate with z G -D(^l) satisfying the differential 
equation (|3.6|). Moreover, z(t) is the unique solution with z(0) = z provided / is such that 

lim [ P \\f(s)\\ds = 0. 



Now, substituting f(s) = F (z(s)) into (|3.12j ) leads to an integral equation 

z (t) = e~ tA z + f e - {t - s)A F {z{s)) ds (3.13) 
Jo 



whose solution, whether it exists, also solves the initial value problem ( |3 . 1| ) provided F (z(s)) is 
shown to be locally Holder continuous on the interval < t < T. 

To formulate the necessary condition on F and state our results, let B 1 = D^A 1 ), 7 > 0, denote 
the Banach space with the graph norm 

ll/ll 7 := \\AV\\. 

F : B 1 — ► L 2 (— 7r, 7r) is said to be locally Lipschtzian if there exist U C B 1 and a finite constant 
L such that 

\\F( Zl ) - F(z 2 )\\ <L\\z l -z 2 \\ 1 (3.14) 

holds for any z^zjE U. 



Theorem 3.2 The initial value problem (\3.dj ) has a unique solution z(t) for all t G IR + with 
z(0) = z G B 1 ' 2 . In 
compact set in B 1 ^ 2 . 



z(0) = z G B x l 2 . In addition, if Wz^t)]]^ is bounded as t —> 00, the trajectories {z(t)} t>0 lie on a 



9 



Proof. The proof of Theorem [O] will be divided into four parts. Firstly, F(z(t)) will be shown 
to be Holder continuous under the Lipschtzian condition ( 3.14] ), which establishes the equivalence 
between the integral equation ( |3.13| ) and the initial problem ( |3.6| ). Secondly, the Banach fixed point 
theorem will be used to show the existence of a unique solution z(t) of (|3.13|) for < t < T. Hence, 
by a compactness argument, the solution z(t) will be extended to all t 6 R + . Finally, assuming 
stays bounded for all t > 0, we conclude the proof. We have to wait till Section || 



that \\z 



1/2 



for the boundedness hypothesis to be established. 

Part I: Continuity. Let us show that F : D(A 1 ^ 2 ) — > LI p (—tc,tc) given by F{z) = —2azz' is 
locally Lipschitz. We note that D(A l l 2 ) = (— it, n) where H% (— tt, it) is the Sobolev space of 
odd periodic functions which have distributional derivatives up to order k. It thus follows that, if 
z G , then z(x) = f£ z'(£) d£ is absolutely continuous with 



SUp \z(x)\ < V 27T ll^ll^ , 

x€[— 7T,7r] 



by the Schwarz inequality. Moreover, using ( [3.10D , we have 

\\F( Zl ) - F(z 2 )\\ < 2a{\\z 1 (z[-z' 2 )\\ + \\(z 1 -z 2 )z' 2 \\} 

< 2aV2n | H^i || H^x - z 2 \\ 1/2 + ||-2i - z 2 \\ ||^2||i/ 2 | 

which satisfies (|3.14| ) with 7 = 1/2 and L = 2a\[2~n ^ || ^1 1| x / 2 + II^Hiy 



(3.15) 



1/2 J" 



Suppose that z : (0,T) — ► £> 1//2 is a continuous solution of (|3.13| ). From the estimate (|3.10| ) 
we have 



(« 



-HA 



^^-Hi/2 < r 

Jo 



< 





Ae~ 


Jo 




f\ 


\A X - 6 


Jo 




Ci-s 


f~ 

Jo « 


Cl - S h s C s 




1 1/2 



A l - s e- sA \\ds \\A s e- TA w\ 



1/2 



(3.16) 



1-5 



ds \\A 5 e- TA w\ 



1/2 



for < 5 < 1/2 which can be used in the equation (|3.13j ) along with (|3.14|) , to get 



\z(t + h)- z(t)\\ 1/2 < || {e~ hA - I) e-^o|| + ft || {e- hA - I) e-^ A F(z(s))\\ ds 



+ jl +h \\e-^ t+h -^ A F{z{s))\\ ds<Kh 5 



(3.17) 



for some constant K < 00 in the open interval (0,T). Combined with ( |3.14| ), this implies the 
Holder continuity of f(t) = F (z(t)) and the equivalence between the equations ( p.6|) and (|3.13j ). 
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Part II: Local existence. Let V = {z G B 1 ^ 2 : \\z — z \\ < e} be an e-neighborhood and let L be 
the Lipschitz constant of F on V. We set B = \\F(zq)\\ and let T be a positive number such that 



(e~ M -')*o|| 1/2 <f (3.18) 



with < h < T and 



C 1/2 (B + Le) [ s' 1 ' 2 e-™ ds < £ - (3.19) 
Jo ^ 



hold. 



Let S denote the set of continuous functions y : [to, to + T] — > B 1 ^ 2 such that \\y(t) — z \\ < e. 
Equipped with the sup-norm 

||y|| T := sup ||y(t)|| 1/2 

t <t<t +T 

S is a complete metric space. 

Defining <%] : [t , t + T] — ► B 1 ' 2 for each y G S by 

r-t 



§[y](t) = e- (t - to)A z + ( e- {t ~ s)A F {y{s)) ds 

Jt 



we now show that, under the conditions (|3.18|) and ( |3.19| ), $ : S — > S is a strict contraction. 
Using 

||iW))ll < \\F(y(t)) - F(z )\\ + ||F(2b)|| < L \\y(t) - z \\ 1/2 + B<Le + B 
and ( |3. ICQ , we have 



my](t) - z \\ 1/2 < \\(e-^ A -l)e- tA z \\ l/2 + WA^e-^W \\F(y(s))\\ ds 

J to 

< - + C l/2 {B + Le) I s~ l/2 e~ cs ds < e 
2 io 



and since is continuous by an estimate analogous to (|3.17| ), G S. 
Analogously, from ( |3.14|) and (|3.19| ), for any y,w G S 



My](t) - *M(t)IU < H^e-^H \\F(y(s)) - F(tt>(*))|| <fa 



to 



< C 1/2 L / s- l ' 2 e- cs ds 
Jo 

holds uniformly in t G [to, to + T] concluding our claim. 



1 M 

\y-w\\ T < - \\y-w\\ T 
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By the contraction mapping theorem, $ has a unique fixed point z in S which is the continuous 
solution of the integral equation (|3.13|) on (t , t + T) and, by Part I, is the solution of (|3.6|) in the 



same interval with z(t ) = z G B^~. 

Part III: Global existence. As the set U where Q3.14| ) holds is compact, the same T can be chosen 
in Part II for any initial condition zq G U. Moreover, if I\ = (t i , ti + T) and I2 = (£2,^2 + T) are 
two intervals containing to, then there exist Zo,i, 20,2 £ U such that the two solutions z\{t) and z%(t) 
of equation (|3.6| ) on Ji with Z\{ti) = z 0> i and on I 2 with 22(^2) — z o,2, respectively, coincide in the 
open interval IiH^- As a consequence, one can define an open maximal interval J max = t + ) 
(containing the origin), where the solution z(t) of (|3.6| ) is uniquely given by patching together the 
solutions Zj(t) on intervals Ij with Zj(tj) = Zqj. By construction, there is no solution to (|3.6|) on 



(to,f) if t' > t + . Therefore, either t+ = 00, or else there exist a sequence {£n} ng N + > with t n — > t + 
as n — > 00 such that z(t n ) tend to the boundary dU of the compact set U . 

It thus follows that, if t + is finite, the solution z(t) blows-up at finite time. In what follows we 
show that ||£(£)||i/ 2 remains finite for all t > to and this implies global existence of z(t) . Let us 
start with the following generalization of the Gronwall inequality. 

Lemma 3.3 (Gronwall) Let £ and 7 be numbers and let 9 and ( be non-negative continuous 
functions defined in a interval I = (0, T) such that £ > ; 7 > and 

C(t)<9(t)+Z [ (t-ry 1 C(r)dr. (3.20) 
Jo 

Then 

C(t)<9(t)+ [ E'{t-T)6{T)dT (3.21) 
Jo 

holds for t G I, where E' = dE 1 /dt, 

00 1 

^W = Er>7TT) (er(7)tT (3 - 22) 

poo 

and T(z) = / t^e^dt is the gamma function. In addition, if 6{t) < K for all t G /, then 
Jo 

C(t) < KE^{t) < K'e^ )T (3.23) 

ZioWs /or some /imie constant K' . 

Proof. If T is an integral operator given by the convolution 

TC(t)=e f (t - rf- 1 C(r) dr , (3.24) 
^0 



then the inequality (|3.20|) can be formally solved by 

00 

at) = 9(t) +J2^ n m 



n=l 
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where T n is also an convolution integral operator which can be explicitly evaluated by the Laplace 
transform, 



(a"(-,))" / l{t-Tre(r)dr = (f n *e)(t) 



r (n7 + l) v v " Jo dt 

with/ n (i) = (er( 7 )r) n /r (n 7 + i). 

Equation ( |3.21|) (and (|3.23|) by the fundamental theorem of calculus) thus follows by setting 
E 7 (t) = J^neN fn(t)- Note that this series is absolutely and uniformly convergent in t G I, with 
E 7 (0) = 1, and it cannot grow faster than exponential 

EJT) ~ -e™ T (3.25) 

7 

as T — > oo (see Lemma 7.1.1 in Q). This concludes the proof of Lemma |3~3 . 

Taking the graph norm of ( [3.13|) , we have in view of (|3.9|) , (|3.10|) and ( |3.25|) 

\\m\i/2 < ||^ ( " o) ^o|| 1/2 + ^ f \\A^e-^ A \\ \\z{s)\\ 1/2 



□ 



to 



< C\\z \\ 1/2 + L (t - s)~ 1/2 \\ z(s)\\ 1/2 ds (3.26) 

J t 

< C exp (LCi/ 2 y/7rt) \\z \\ 1/2 , 

which is finite for any t G 1R+. 

Part IV: Compact trajectories. Since B 1 C B 1 ^ 2 has compact inclusion ifl/2<7<lQ, it 
suffices to show that ||2;(i)|| remains bounded as t — > oo. The hypothesis ||^(i)|li/2 < 00 combined 
with ( [3.15| ) implies the existence of C < oo such that, analogously as in Q3.26j ), 

\\z{t)\\ 7 < \\e- tA z \\ 7 + £ H^e-^H \\F(z(s))\\ ds 

< C 1 - 1/2 t 1,2 -~ < e- ct \\z4 l/2 + C' [ {t - s)' 1 e~ c{t ~ s) ds , 

Jo 

which is bounded for t > provided c > (i.e. infA o~(A) > ). Although the spectrum of A is not 
positive if (3 < 8n, we shall see in Section |S| that A in the integral equation ( |3.13| ) can be replaced 
by a positive linear operator L (see Theorems |5.2| and |5.3| ). 



This concludes the proof of Theorem 3.2 



□ 



It follows by analogous procedure that if z\ and z 2 are solutions of ( |3.6| ) differing by their initial 
value in B 1 / 2 , then 
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M*)-* 2 (0lli/ 2 < ||e^ A Ki-^o,2)|| 1/2 + / WA^e-^W \\F( Zl (s)) - F(z 2 (s))\\ ds 



< \\e- tA (z 



i-z , 2 )\\ l/2 + C l/2 L [\t - S y 1/2 e- cs ds \\ Zl (s) - z 2 (s)\\ 1/2 
Jo 



which implies, by the Gronwall inequality, the continuous dependence of z(t) with respect to its 
initial condition. 

We may also consider the dependence of z with respect to the parameter a = /3/(4ir). The next 
statement is a corollary of the above analysis. 



Theorem 3.4 The solution z(t) : M + x B 1 ! 2 — > B 1 ! 2 to the initial value problem ( $.dj ) as a 
function of the bifurcation parameter a and the initial value z is continuous. 

Remark 3.5 It can be shown (see JT^/j that for any initial value zq e B~> ' , < 7 < 1, the solution 
is actually in D(A) at any later time. Moreover, since F : £> 1//2 — > L 2 p (— n, it) is C°° (has Frechet 
derivatives of all orders), it can also be shown that (a, z ) E W + x B 1 ! 2 — ► z(t; a, z ) is C°° for all 
t > 0. 



Remark 3.6 Under minor modifications, one can show existence, uniqueness and continuous de- 
pendence of (\3. 1\ ) in Sobolev space H^ p (—n , n) with norm 1 1 ^ 1 1 x = ||^'|| X 2 (just include the linear 
term of (3~1) in the definition of F). The same results hold for equation in the Sobolev space 

of even and periodic function H^ p (— 7r, tt) with both norms || - 1 1 x and || • || x / 2 - Note from item 2. after 
([3.Q) and ( \3. that a \\z\\ l = H^H^ + 2 H^H^a so, both norms are equivalent. 



4 Equilibrium Solutions 

Time independent (equilibrium) solutions of ( |3.1|) are odd solutions of the ordinary differential 
equation 

a(^"-2^') + 2 ^ = 0, (4.1) 

with periodic conditions ip(—7r) = ip(7i) and n) = "^'{it), a = (3/ (Ait) > , which can be 
written as 

w' = 2p{w — a~ 1 ) 

(4.2) 

p' = w , 

by setting p = ip and w = ip'. 

In this section we give a qualitative and quantitative description of the solutions of (|4.2|) in the 
phase space M 2 and study their implications for the equilibrium solutions of (|3.1|) . Our results are 
summarized as follows. 
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Theorem 4.1 The stationary equation (4-1 ) has two distinct regimes separated by a = 2 (f3 = 87r). 
For a > 2, ip = is the unique solution. For a < 2 such that 2/ [k + l) 2 < a < 2/k 2 holds for 
some k G N + , there exist 2k non-trivial solutions i[)j~,ifjj , j = 1, . . . , k, with fundamental period 
2n/j, ipj (—x) = —ipj(x) and ip7(x) = ip^x + tc). Moreover, each pair of non-trivial solutions 
bifurcate from the trivial solution ipo at aj = 2/j 2 (f3j = 8tt / 'j 2 ) with lim ip, = 0. 

In the phase space, these solutions (^j'^ipj), are closed orbits around (0,0) whose distance from 
the origin increases monotonically as a decreases. Numerical computations indicate that these 
orbits approach rapidly to the open orbit {(a -1 , a -1 x) , x e M.} from the left as a — » 0. 



Let us begin by stating the general properties derived by the same tools used in the analysis 
performed in Section [3]. 

The vector field / : M 2 — ► M 2 , 

(w,p) — ► f(w,p) = (2p(w - a' 1 ), w) , 

in the right hand side of ( f4.2| ), defines a smooth autonomous dynamical system. It thus follows 
from Piccard's theorem (see e.g. [ |CL|| ) that there exist a unique solution (w(x),p(x)) of this system, 
globally defined in M. 2 , with (w(0),p(0)) = (w ,p ). As we have seen in Section |3|, the existence 
of a global solution and its continuous dependence on the value (wq,Po), and on the parameter a, 
follow from Gronwall's lemma, which holds here in its standard form. As a consequence, the phase 
space IR 2 is foliated by non-overlapping orbits 

lp = {(w(x), h{x)) : x G R and P = («;(0),p(0))} 

which passes by P = (w ,p ) G M. 2 at x = 0. Note that, by varying continuously P and a, the orbit 
7p varies continuously in the phase space. 

We shall now determine the values (P, a) by which the solution of ( |L2| ) defines closed orbits. 
Note that the orbits are symmetric with respect to the u>-axis, L = {(w,0) : w G IR}, since the 
system of equations ( |4.2f ) remains invariant if the sign of both, x and p, are reversed. As we shall 
see, there is no loss of generality if the initial value (w(0),p(0)) = P belongs to L. We write 

IP = Iwo- 



Proposition 4.2 Every orbit 7p is determined by a single value P in the positive semi-axis L + = 
{(w ,0) : w > 0}. For w > 0, the orbit j Wo is either closed or unbounded depending on whether 
auiQ < 1 or auiQ > 1, respectively. The orbit j a -i = {(a~ 1 ,a~ 1 x) : x G M} separates the phase 
space R 2 in such way that jp is closed if P is on the left of 7 a -i and unbounded otherwise. In 
addition, if w = 0, then 70 = {(0,0)}, and the origin is enclosed by every closed orbit. 



Proof. The proof of Proposition follows from an explicit computation. By the chain rule, 
equation (f4.2|) can be written as 



dp 
dw 



w 



2p (w 



or 



(4.3) 
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provided aw ^ 1. The trajectories j m , obtained by integrating 2pdp = wdw/ (w — a x ) with 
initial point P = (w , 0), 



are portrayed in Figure 1. 



2 , / 1 - aw \ 

p = w — w + a In J 

VI - atwoj 



(4.4) 




w 



Figure 1: Trajectories of the dynamical system 
We note that P = (0, 0) is the only critical point of fl4.2|) which is a center for all a > since, 



by linearizing f(w,p) around P = (0, 0) gives a matrix whose eigenvalues are A± = ±iV2ar I . This 



implies that 7 = {(0,0)} and the orbits r f m with w sufficiently closed to are, in view of (|4.4j ) 
ellipses defined by the equation 2a~ 1 p 2 + w 2 = C. 

When aw = 1, using mathematical induction and equations ( |4.2|) with (u;(0),p(0)) = (w ,0) 
we have 



dx n 



(0) = 



for all n > 1, which leads 



j a -i = {(a , a x x) : x G M} . 
Hence, if u = oo(P) denotes the set of limit points (the u -limit set) given by 
co(p) = \(w*,h*) e R 2 : lim (w(x n ),h{x n )) = (w*,h*) 



(4.5) 
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for some sequence of points {x n } such that x n — > oo , as n — > oo, j a -i separates two different type 
of orbits: to{P) = 7p or u(P) = {00} depending on whether the point P is at the left or at the 
right of 7 a - 1. 

□ 

Proof of Theorem |4.1| . The stationary solutions satisfy (|4.2|) with periodic conditions w(0) = 
w(2ir) and p(0) = p(2n). By fixing the period T of an orbit ^y wo in 2tt, the label Wq becomes 
implicitly dependent on the parameter a. In view of Proposition [4.2| , Theorem |4.1| follows if for 
a > 2, except by the orbit 70 = {(0,0)}, no (non-trivial) solution has period T = 2n and for 
a < 2 there is a one-to-one correspondence between wq and a for T fixed at any value 27r/k, 

fc = l,... 



More precisely, let T = T(a, wq) denote the period of the dynamical system (|4.2|) with initial 
value (w(0),p(0)) = (w ,0): 

T= f dx = 2 / ^ , (4.6) 

K J w 

where, by symmetry, the second integration is over the semi-orbit above the w-axis. For T> = 
{(a, Wo) G M + x M + : aw < 1}, we set 

J 

and note that Gj : T> — > K is a continuous function of both variables satisfying 

Gj {2/f, 0) = 0. (4.7) 

To see Q4.7D , we compute the period T L of an elliptic orbit, e.g. {(2/a)p 2 + w 2 = 1}, of Q4.2D 
linearized at the origin (f(w,p) replaced by (2a _1 P) w ))) 

T L if^- £, -f, (4.8) 

h VI - (2/a) P 2 V 2 ' 

and note that lim„, ^o ^(«, ^o) = ^l- Continuity follows from the general properties stated previ- 
ously. 

Hence, provided 

dT 



dw ( 



> (4.9) 



holds for all (a, wq) G by the implicit function theorem, there exists a unique (strictly) mono- 
tone decreasing function Wj : [0, 2/j 2 ] — ► R + with Wj(2/j 2 ) = such that Gj(a,Wj(a)) = 0. Note 
that (fO|) and 

T(a,wo) = \/aT(l,awo) (4-10) 

imply that T is an increasing function of both a and wq, independently. This fact, which can be 
seen by rescaling fl4.2| ) by x — > x = x/ \/a, w — > tZJ = aw and p — > p = v^P; explains the monotone 
behavior of Wj. 
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It thus follows that, if a < 2, for each j = 1, . . . , k such that 2/ (k + l) 2 < a < 2/k 2 holds, a 
unique function Wj such that Wj(2/j 2 ) = exists. The non-trivial solutions ip^ 1 , . . . , ip^ of ( (4.1| ) are 
the p-component of 7^., j = 1, . . . , k, which winds around the origin j-times: is 27r-periodic 

odd function with fundamental period 2ir/j, (^//(O) > and satisfies ip^{x + 7r) = ijjj{x). If 
a > 2, because T(a,w ) is a strictly increasing function of w and T(a,0) > 2-n (see eq. ( ^8|) ), 
there is no solution of Gj(a,w ) = besides Wj(a) = for j = 1. This reduces the proof of 
Theorem JO] to the proof of inequality ( f4.9|) . 

To prove (|4.9|), it is convenient to change variables. Let 



q — In (1 — aw) 



(4.11] 



be defined for aw < 1. From ( |4.1U| ), there is no loss of generality in taking a = 1. The system of 
equations (|4.2j ) under this condition is thus equivalent to the following Hamiltonian systemP] 



q' 



p 



2p 



(4.12) 



whose energy function is given by 



The trajectory equation 
level equation H (q, p) = E with 



H(q,p) =p 2 + e q -q-l. (4.13) 
when written in terms of the g-variable, gives exactly the energy 



E = -w - In (1 -w ) . (4.14) 
We denote by 7^ the orbits of ( f4.12| ) and note that, in view of the fact 

dE = w 
dwo 1 — wo 

there is a one-to-one correspondence between the two families of closed orbits {7 W0 , < w < 1} 
and {-ye, < E < 00}. 

Now, let T = T(E) be the period of an orbit je, 

" q+ dq 



T 



dx 



IE 



P 



Using the energy conservation law, we have 

p = p(q,E) = 
where the potential energy is given by 

v(q) = e q 



E - v(q) 



q-l 



and q± = q±(E) are the positive and negative roots of equation v (q) = E. 

df 



Equation ( |4.9| ) holds if and only if - — > holds uniformly rn.EE 

dE 

the monotonicity criterion given by C. Chicone [0 (see also [ |CG|| ' 



(4.15) 

(4.16) 
(4.17) 

But this follows from 



We thank G. Benfatto for explaining this tranformation and for pointing us equation (iA) in a footnote of 
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Lemma 4.3 Let v G C 3 (R) be a three-times differentiable function and let f(q) 
force acting at q. Ifv/f 2 is a convex function with 



v 
P 



6v (v"f - 3 (v'f v" - 2vv'v'' 



> 



-v'(q) be the 



(4.18) 



then the period T is a monotone (strictly) increasing function of E. 
Proof. It follows from (|4.16|) two basic facts: 



dp = f_ 
dq 2p 



and 



p(q±,E) = 0. 



(4.19) 



These will be used for deriving an appropriated integral representation of dT/dE. 
Let 



/ rip) <!<1- 



Integrating twice by parts, gives 



3 V/ 2 



9+ 



1+ 



pv 
2? 

r 



+ J (Pf)' J 2 dq 



— + vp^pr ) dq 



2j q _ \2p 



(4.20) 



in view of (|4.19|) . Note that f(q±) ^ since 

v'(q±) = v(q±) -q± = E-q± 

vanishes only at E = 0. This follows from the fact that v is a convex positive function with v(0) = 
and asymptotic behavior v(q) ~ q — 1 and ~ e aq , as q goes to — oo and oo. 

Now, using iy/f)' = v' I f — vf'/f 2 = —1 — vf'/f 2 , and integrating by parts, we continue 



K 



'<?- 
1 ri+ 



Yp- p \-f) 'Pi dq 



l- p ) dq ~\p (j 



(4.21) 



where in the last equation we have used v = E — p 2 . 
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From (gjjD , (ggg) and flOID , we have 



ET = 2 pdq + 



1+ 



V 



3.L r [T 2] dq 



Differentiating this with respect to E and using ( 4.19| ), gives 



~ f q+ dq f q+ 



V 

, V J,, \f- 



dq 



which, in view of ( |4.15|) and the assumption of Lemma [4.18| , implies 

df 1 ru 



dE E 



P {jj J (hl - > ■ 



□ 



It remains to verify ( (4.18| ) for v given by ( |4.17| ). By an explicit computation (see Chicone [|Cj]) 



js) K) 4 = <■".'/<</> 



where 



g (q) ■= e 2q + 4 (1 - q) e q - 2q - 5 



is such that g(0) = g'(0) = and g"(q) = Ae q v(q) > 0. This implies g(q) > (g(q) = only if 
q — 0), the hypothesis of Lemma |4.3| and concludes the proof of Theorem |4.1| . 



□ 



Turning back to the Coulomb system problem, some remarks are now in order. 



Remark 4.4 Recalling v(t,x) = u x (t,x) and denoting X* = lim A n the charge activity at the fixed 
point, we have from ( \2.1b]) 



q& I q&L 



and 



^'(0) = ]Tg 2 A*(g) /$>*(g) >0 



q& 



These boundary conditions select ip~^ , j = 1, . . . , k, as being the only physically meaningful station- 

ary solutions and implies 4> + {x) = / ?p + (y) dy > on (— 7r,7r). 

Jo 
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Remark 4.5 The value a = 2 is a bifurcation point as one can see by linearizing d4 ■ A ) about 
■0 = 0. The linear operator L[0] = A given by (\3. TO in the subspace of odd 2ir -periodic functions 
has eigenvalues and associate eigenf unctions given by ( \3.8{ ). Hence, if a > 2, the eigenvalues are 
all positive and ip = is locally stable. When a < 2 (but close to 2) a single eigenvalue becomes 
negative and one can apply the Crandall-Rabinowitz bifurcation theory /(£]/ to locally describe the 
stable solution which bifurcates from the trivial one. Note that Crandall-Rabinowitz theory can 
also be applied in the neighborhood of aj = 2/j 2 , j > 1, in the orthogonal complement of the 
span sinmrr, m — 1, — l} corresponding to the odd functions with fundamental period 

T = 2tx / j . These points were referred to in the introduction as a sequence of instability thresholds. 



In Theorem LI we have given a global characterization of the non-trivial stationary solutions. 



Remark 4.6 In the sine-Gordon representation, the effective potential (p(x) = ip(y) dy = 
x 2 j (2a) at 7 a -i corresponds the Debye-Hiickel regime with Debye length a. Although this regime is 
not reached for all (3 > 0, it gets closed quite fast as (3 = Ana approaches 0. Numerical calculation 
is shown in Figure 2. Note that at a = 1 ((3 — 4ir), Wi cannot be distinguished from a' 1 (numerical 
error is in the sixth decimal order). 




—• ■ ■ ■ ■ ■ ■ — 1 ■ ■ ■ 1 ■ ■ ■ 1 ■ ■ ■ 1 ■ ■ ■ 1 — a 

0.8 1.2 1.4 1.6 1.8 2 



Figure 2: Comparison between the initial value function for the periodic orbit of period 2ir, wi = 
wi(a), and for the nonperiodic Debye-Hiickel orbit, WDH(a) = a" 1 - 



Remark 4.7 The derivative of ( K.dj ) with respect to Wq, computed from equation ( \4-4 ), 



dT 



2awc 



dw 1 — aw 



sign [w) 



2 (1 — aw)p 2 /a 
[1 + 2(1- aw) p 2 /a) 



dp. 



indicates that an estimate from below can be very delicate to obtain. Note sign (w) changes along 
the orbit r y Wo . This shows how amusing Chicone's monotonicity result is for the problem at hand. 
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5 Stability 



Let z(t; zq) denote the solution of the initial value problem (|3.6|) - (|3.7|) . It follows from the analysis 
in Section [3] that 

S(t)z = z{t;z ) (5.1) 

defines a dynamical system on a closed subset V C D (A) of B 1 ^ 2 with the topology induced by the 
graph norm IHI^- Note that z(t; z ) is continuous in both t and z with ^(0; z ) = z and satisfies 
the (nonlinear) semi-group property S(t + r)z = z(t; z(r; z )) = S(t)S(T)z . 

This section is devoted to the stability analysis of the equilibrium solutions described in Section 
|4]. By local stability it is meant that z(t; z ) is uniformly continuous in V for all t > 0: given e > 0, 
\\z(t; Zq) — z(t; -Si) || !/ 2 < e f° r alH > and Z\ G V such \z\ — ^o|li/2 < ^ f° r some 5 = 5(e) > 0. It 
is uniformly asymptotically stable if, in addition, lim ||2:(t; z ) — z(t; zi)\\ 1/2 = 0. 



t—tOO 



The Liapunov (global) stability analysis as developed by LaSalle and applied to semilinear 
parabolic differential equations by Chafee and Infante ||C]|| (see also |H[]) will also be discussed and 
extended in this section. 

Let us begin with the local analysis. 

Theorem 5.1 (Local Stability) There exists a neighborhood Li C £> 1//2 of the origin such that, 
if a > 2 and z in U , then ip = is stable, i.e., lim \\z(t; ^o) II 1/2 = 0. If a < 2 is such that 

2/ (k + l) 2 < a < 2/k 2 holds, among all equilibrium solutions of ipo, ipf, j = 1, . . . , k, ipf 

are the only asymptotically stables. So, there exists p > such that if \\z — "01 1 1/2 < P, then 
lim \\z(t; zq) — ^||i/ 2 = for vp = ipf and, for any sequence {z n } >l with lim \\z n — ip\\ = 0, we 

1/2 > e > f or a U n an d ip = ipj 



have sup \\z(t; z n ) — ip\\ 1/2 > e > for all n and ip = ipf, j 7^ 1 



It is convenient to consider the equation 



for ( = z — ip where ip is a solution of (|4.1|) . Here 



L( = L [ip] C = -otC" + 2aip(' - 2 (1 - o^') C (5.3) 

is the linearization of the differential operator (ft.l|) around ip an d F is as in ( |3.7| ). Note L = A and 
(U) reduces to if = ^0 = 0. 

Proof. The proof of the Theorem |5.1| follows from the next two theorems. 



Theorem 5.2 If the spectrum o~(L) of ( \5. 3j) lies in {A G K. : A > c} /or some c > ; i/ien £ = 
zs i/ie unique uniformly asymptotically stable solution of ( \5.2[ ). On the other hand, if o~(L) H 
{A G R : A < 0} ^ 0, &en C = unstable. 



22 



Theorem 5.3 Let L = L[ip\ be given by (\5.9{) . Then cr(L) > whenever if> = ip Q and a > 2 
or ip = ipf and a < 2. If a is such that 2/ (k + l) 2 < a < 2/k 2 holds for some k e N + , then 
a(L) n{Ael:A<O}^0 for ijj = ip Q and $ = 1$, j = 2, . . . , k. 



Proof of Theorem |5.2| . We shall prove only the first part of Theorem [|72] and refer to Theorem 
5.1.3 of Henry's book |H| for the instability part. 

It follows from ( |3.13| ), ( |3.1U| ), ( |3.15| ) and the hypothesis on o~(L) that 



||CWII 1/2 <^i/ 2 e- ct ||Co|| 1/2 + e / (t-s) 



-1/2 -c(t-s) 



IIC(s)|| 2 /2 ds 



(5.4) 



with c> 0, Ci/2 = l/\/2e and £ = 2v / 2~7ra. 

Let us assume that ||C( s )lli/2 — P on a interval (0,t) for some p satisfying 

If || Co II 1/2 — P\j\ 1 then equation (|5.4|) can be bounded as 

nc(t)ii 1/2 <f + P 2 e^ (*-*r 1/2 c- 



(5.5) 



i. e., p < 



4ira 



- ^-V2 „-c(t-«) 



(5.6) 



and this implies the existence of a unique solution of ( ^2|) with ||C(*)|li/ 2 — P f° r a ^ £ > 0. Note 
that || Co II 1/2 < P an d ^ ^i i s the maximum value under which ||C(^)ll i/2 < P f° r all < t < ti, then 
either ||C(^i)lli/2 = P or *i = 00 • But the first case is impossible by ( |5.6|) . 

Going back to (|5.4j), using || C( s )lli/2 < P an d a slightly modification of Gronwall inequality 



]f) with £i /2 (t) = {p£Vnt 1/2 ) n /r(n/2 + 1) , we have 



n=0 



IICWIIl/2 < Ci/2||Co|li/2^1/2(t) 



-ci 



< 



(VilCoiU (i + p^/ 2 ) e -( c -^)* 



< 2illCo|| 1/2 (l + 



1 /ct 



,-3ct/4 



2 V 7T 



in view of fl5.5|) . This proves the stability statement of Theorem |5.2| , since (|5.2|) defines a dynam- 
ical system in a closed subset V p = |c € -B 1 / 2 : HCII1/2 — p} "with hm ||C(X)Hi/2 = if 1 1 Co 1 1 1/2 = 



F0 



11/2 <PV 2" 



□ 
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Remark 5.4 One can actually show that if c = inf^cr(L) then ((t; Co) = z(t; Zq) — if> decays expo- 
nentially fast to as 

C(*;Co) = ^(Co)e- ct + £(t;Co) 

where ||e(£; Co) II 1/2 — @ IIC0II1/2 e ~ c '* w ^h < c < d and k : V p — > N (L — cl) is continuous 
and such that k{0) = 0, where N (L — cl) is the one-dimensional span of the eigenfunction of L 
associated to c. 



Proof of Theorem |5.3| . Since L[if)o] = A, Theorem |5.3| for if) = ip with a > follows from the 



spectral computation in ( |3.8| ). 

Now, let if) be a nontrivial solution of the equilibrium equation ( f4.1|) and note that if>(0) = 
if) (it) = by parity. According to Theorem |5.2| , if> is asymptotically stable if o~(L) > and unstable 
if a(L) n {A < 0} ^ 0. 

Let ip be the solution of 

L[i/f]<p = (5.7) 

in the domain < x < tt satisfying 

y>(0) = and <p' (0) = 1. (5.8) 

As in JJ], we shall use the comparison theorem to establish that if) is asymptotically stable if 
(p(x) >0on0<x<7r and unstable if <f(x) < somewhere in < x < tt. 



To apply the comparison theorem and complete the proof of Theorem |5.3|, let 

p( x ) :=e -2/ ^(v)di/ (5.9) 

be the weight which makes L a self-adjoint operator: 

p L[iP}( = -a (pC) 1 - 2p (1 - aif)') ( . (5.10) 

Note that (L(,r]) p = ((,Lr]) p for any odd periodic functions ( and rj of period 2tt were (f,g) p '■ = 
f(x) g{x) p(x) dx. 



Theorem 5.5 (Comparison) Suppose £i and (2 ore two real solutions on the domain (0, n) of 

pL[ij;]C = fi, z = l,2, 

respectively, with Ci(0) = Ci(vr) = 0, Ci(0) > and ( 2 {0) = 0, C 2 (0) > 0. If (i > and = fi{(;x) 
is such that 

h > fi (5.H) 
on (0, 71"), i/ien (2 must vanish at some point of this domain. 
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Proof. Let assume that (2 > on (0, 7r). Then, from ([5.10 ) and the hypotheses of Theorem [5]5 
we have 

2 f (/ 2 - /1) ^ = (Ci, LC 2 ) p - (LCi, C 2 ) p 
./o 

= 2a [* [(p( , 1 )'C2-(i(pQ'} dx 



= 2a ^ [p (([(2 - Ci Q}' dx 
Jo 

which, in view of the boundary conditions and (|5.11 ), implies a contradiction 



>0. 

Note that ([(n) < since £i > on (0, ir) and Ci( n ) — 0- So, there must exist x G (0, 7r) such that 

= 0. 

□ 

If we consider the eigenvalue equation 

L[ij>]9 = X9 (5.12) 
on (0, 7r) for the smallest eigenvalue A in the space of odd periodic function, 9 satisfies the conditions 



of (1 in Theorem ^75] with fi = Xp(. Note the eigenfunction associated to the smallest eigenvalue 
may be chosen to be positive in the domain (0, n). 

Applying Theorem |5.5| for ( |5.7| ) and ( |5.12| ) we arrive to the following stability criterium: 



Criterium 5.6 The smallest eigenvalue A of L [ip] is positive if if > on (0, 7r) and negative if 
there exist x G (0, 7r) such that (f(x) = 0, where ip is the solution of equations ( \5. 7j ) and ( \5.8j ). 

Now, for a given non-trivial stationary solution ip let 

X = c(-ar + ^) , (5.13) 

where c > is chosen so that = 1- ^ follows from the equation —aip" = 2(1 — aip') ip (see 
LTD), that 



X (0) = and x > 

whenever ip > (recall ^(0) = and 1 — aip' > for all closed orbits). Moreover, we have 

Proposition 5.7 

L[ij}] x = 8co 2 ^ (ip') 2 > (5.14) 
on the same domain (0,x) that ip > 0. 
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Proof. Differentiating (|4.1|) twice, 

-a {if)")" = ~2aip + 2(1- 3aif)') ^" , 

and using ( |4.1| ) again, gives 

L[V#" = -a(^")" + 2a^(^") / -2(l-a^ / )^" 
= -iocif/tf/' 
= 8 (1 - an//) if) if)' 



In addition, we have 



-oaf/' + 2aipip' - 2 (1 - a^') ^ 
2aif)if)' 



which combined with the above equation, gives the equality in Proposition 5.7. 



□ 



Completion of the proof of Theorem |5.3| . We are in position to prove Theorem |5l| for non- 
trivial equilibrium solutions. Let x De given by ( |5.13| ) with if) = if)f. Then x > on (0, 7r) and 
Theorem |0| can be used to compare equation (|5.14j) with G5.7|) . This yields <p > x > on (0, it] 



which implies the stability of if>i by Criterium [5]6 
For instability, we observe that if)' satisfies 

L[ifj}ip' = -aif)'" + 2aif)if)" - 2 (1 - aif)') if)' 

= (-aif)" + 2aipip' - 2if))' = , 

in view of equation (|4.1| ). Recall that if) = if)^ with j > 2, has fundamental period 27r/j and satisfies 
I j) — '0 / '( 7r / j) = by the odd parity and equation (JO). Since if)'(0) > 0, this implies if) < 

on (ir/j, 2iv/j) and the minimum of if) is attained at x — — r- Since if)' and <p satisfies the same 

2j 

self-adjoint equation pL[if)]( = 0, their Wronskian 



W^,if)';x) 



-apif' 



ifj' 
-apifj" 



= ap ((p'ip' - (pip") = aif)'(0) > 

is a non-vanishing constant (recall p(0) = 1, y?(0) = and ("0/) (0) > 0). As a consequence 

W (ip, if)'; ir/j) = -ap(x)<p (x) if)" (x) > 

implies ip(x) < because ip" (x) > 0. It thus follows from the stability criterium that ipj~, j = 
2, . . . , k, are unstable since x G (0, n) provided j > 2 and there exist x G (0, 7r), x < x, such that 
ip(x) = 0. 
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By a slight modification of these arguments, one may conclude the stability of ip l and instability 
of ip~ , j ' — 2, . . . , k, as well. This concludes the proof of Theorem |^ and, consequently, the proof 



of Theorem FTl 



□ 

Now, we turn to the Liapunov stability analysis with a proof of global stability of the trivial 
solution O = 0. 

Let V be a real-valued functional on the subspace of absolutely continuous function of D(A) 
given by 

{(V 1 - v')\n(l-av')+v' -v 2 } dx (5.15) 

•TT 

and notice that V(0) = and V (77) = W(rj) + o (||^|| 2 ), as 1 1 77 1 1 — > 0, where 

W{v) = \ f (av' 2 -2v 2 ) dx, 

by Taylor expanding g(w) = (a -1 — w) In (1 — aw) — w around w = 0. Observe that W (v) = 
(1/2) |Hli/ 2 if a > 2 and since g{w) — (a/2)w 2 > if aw < 1, V (v) < W (v) holds on the space 

V={veHl p nHl p :av' < 1} , 

of odd, positive and 27r-periodic functions with distributional derivative up to second order. 
A Liapunov function V of a dynamical system {S(t),t > 0} satisfies 

V(v) = hmi (V (S(t)v) - V(v)) < (5.16) 



for all dgV. We now show that (|5.16j) holds if S(t) is given by equation (|3.1|). More precisely. 



Proposition 5.8 Let S(t)vo = v(t; vq) be the dynamical system in V given by (V5. Then, the 
pair of functions 

p{w) = (5.17) 
1 — aw 

and 

$(p, w) = (a -1 — w ) In (1 — aw) + w — p 2 (5.18) 
generate the Liapunov function given by ( \5.15 ): 



V(v)= / ^{v,v x )dx with V(v)= / p(v x )v 2 dx . (5.19) 
Jo Jo 
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Proof. Note that, from the parity of v the integral in ( |5.15 ) can be made over [0, 7r]. By the 
calculus of variations and equations ( |5.17| ) and ( |5.18| ) we have 



V(v(t,-)) 



d <9$ <9$\ , <9$ 
^-)v t dx+ -—v t 

ydx OV x OV J OV x 



[ ( — — ln(l — av x ) + 2v) v t dx 
Jo \ dx J 



p{y x ) (av xx — 2avv x + 2v) v t dx . 



(5.20) 



where v t (t, 0) = v t (t,n) — 0, t > 0, in view of the boundary conditions on V. Since p{w) > for 
aw < 1, this with (|3.1| ) concludes the proof of Proposition. 



□ 



Remark 5.9 We have used the construction method based in the Euler-Lagrange equation to find 
this Liapunov function (see e.g. Chap. 2 of Zelenyak, Lavrentiev and Vishnevskii \ZL\\j ). A 
sufficient condition for ^5.1^\ ) hold leads to a first order partial differential equation for p 



wp p (1 - aw)pp w = -2pp 

a 

whose characteristics are given by the orbits 7 WQ described in Section^ in the study of the equilibrium 
solutions of ( \3. 1\ ). Note that equation ( \5.18j ) is the Lagrangian associated with the Hamiltonian 
( j.lS^ j (with q defined by ( j.ll )). Due to the requirement aw < 1, our particular solution takes 
into account only the closed orbits. There may be other suitable choices which includes all orbits. 



The proof of global stability of 0o requires that a subspace of V be invariant under the flow 
equation ( |3.1| ). This is shown in the following by using the maximum principle. 



Theorem 5.10 If v(t,x) is a classical solution of equation with initial condition t>(0, x) 

v (x) G V, then av x (t,x) < 1 and 



a 1 (x — 7r) < v(t, x) < a 1 x . 



(5.21) 



hold for all t > and < x < n. 



Proof. Denoting 



L[v] := F(v xx ,v x ,v) - v t , 



(5.22) 



where F(ai, 02, (23) — a(a\ — 2a2 a^) + 2a^ is a continuous and differentiable function of its variables, 
the differential equation (|3.1| ) can be written as 



L\v] = 0. 



(5.23) 



28 



For v satisfying ( |5.23|) with v(t, 0) = v(t, ir) = , < t < r, and initial data v (0, •) = v , let us 
suppose z = z(t,x) and Z = Z(t,x) are such that 

L[Z) < < L[z) (5.24) 

for all (t, x) in D = (0, r) x (0, vr) with y) < < Z(t, y), y = 0, ir and < t < r , and 

^(0, x) < vo(x) < Z(0, x) , 

for < x < 7r. Then, by the maximum principle (see [PW|| , Theorem 12 in Chap. 3), 

z{t, x) < v(t, x) < Z{t, x) (5.25) 

m~D= [0,tt] x [0,r]. 

The lower limit function z is given by 

=9{x-ir) , (5.26) 

with > 0. From (|5\23| ), 

= 20(a0-l)(7T-ar)}, 

is always positive provided ad > 1. 

Analogously, the upper limit function Z is given by 

Z(x)=6x, (5.27) 

from which 

L[Z] = -25(a5 - l)x 

is always negative provided aS > 1. 

Since (|5.24|) holds uniformly in r, equation ( 5.25 ) holds for all (t, x) in R + x [0, ir]. Note that f 
remains bounded irrespective of av x < 1. However, if this condition holds for t — 0, it remains for 
all t > 0. To see this, observe from the equation t> f = v xx + (1 — av x )v with v xx = 0, that the rate 
by which \v\ increases tends to zero when the inequality saturates. The inclusion of the Laplacian 
only smooths v and prevents, even more, v x to increase beyond the threshold. The same argument 
justify the strict inequality (|5.21|) . 



This concludes proof of Theorem 5.10 



□ 



We pause to discuss some properties of the classical solutions of equations (|1.1|) and (|3.3|). Recall 



that u(t,x) = / v(t,y)dy with v satisfying (|3.1|) . 
Jo 

Remark 5.11 Note that the cone C = [u 6 PI : u > 0, au xx < 1} zs invariant under the 
unnormalized evolution For this, let 



M[u] := a(u xx — u 2 x ) + 2u. 
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Ifu(t,x) is a classical solution of (\1.1\) with initial value uq G C, since M[u] = foru = ; we have 
by Theorem 7 in Chap. 3 of JPTW (see also Remark (ii) after this) that u(t,x) > for all t > 0. 
This, however, does not imply that u(t,x) remains positive (recall (\3.5]) ). A proof of this assertion 
goes as follows. 

Theorem 5.1Q implies u(t,x) remains bounded, and u xx (t,0) bounded from above, ifu satisfies 
^3.3\ ) with initial condition uq satisfying a~ l [(x — vr) 2 /2 — vr 2 /2] < uq < a~ l x 2 /2 (by integrating 
( \5.21\ )). The comparison principle applied directly to equation ( \3.3\ ) leads to ( \5. 24 ) with L replaced 
by M and replaced by au xx (t,0). An upper and lower solutions, z and Z, can be obtained from 
the solution of the equilibrium initial value problem (\j.2\ ): 

<$> ± (a,w ;x)= ^ ± (a,w ;y)dy, 



where \I/ ± (q;, wq; x) is the p- component of the closed orbit •y± WQ starting at (±u> ,0). Note that $ ± 
is an even periodic function of period T = T(a,Wo) given by a monotonically increasing function 
of both wq and a with T — > oo as aw f 1 for $ + and as w — > oo for $~ (see proof of Theorem {j.1\ 
for details). $ +( ~) is also a monotone increasing (decreasing) function of x in [0,T/2] and satisfies 

M[$ ± ] = ±aw , 

with $ ± (0) = ^'(O) = $ ±/ (T/2) = and $ (0) = Wq. The lower limit function z is given by 

z(x) = 9 ($~(a, w ; x + x) + ^w ) , (5.28) 

where 9 < 1, w > w , x = x(a,w ,w ) is such that z(0) = 0, with w and w so that T is very 
large and z'(tt) = which can always be done in view of the properties o/$ _ . The upper limit 
function Z can be written also as (\5.2fy ) with $~ replaced by $ + and the second term with minus 
sign. We have 

M\W±] = a9{(l - 9){$ ± 'Y t (w - w )} , 

with W + = Z and W~ = z. In order inequality ( \5.29j) holds uniformly in t, u xx (t, 0) has to remain 
bounded from above and below. Since u xx (t, 0) < a -1 by Theorem \5.1Cj , one may choose 9 arbitrarily 
small in ( \5. 28[ ) and take Wq and wq so large that 9(wo — Wq) > a -1 . In the limit as 9 — > we have 
u(t,x) > and < u xx (t, 0) < a -1 for < t < r, uniformly in r, implying u(t,x) > for all 
t > 0. 



LaSalle's invariance principle allows us to apply Liapunov function techniques under milder 
assumptions. A subset K, C V of a complete metric space V is said to be invariant (positive 
invariant) if, for any Vq G /C, there exist a continuous curve v : K — > K, with v(0) = Vq and 
S(t)v(r) — v(t + t) for allt > and r G R The following two theorems express the content 

of this principle. 

Theorem 5.12 Suppose vq G V is such that the orbit 7(^0) = {S(t)vo,t > 0} through vq lies in a 
compact set in V and let u(v ) denote its u -limit set, i.e., 

T>0 

(see ( \4-o] ) for alternative definition) . Then uj(vq) is nonempty, compact, invariant, connected and 
dist (S(t)u , u(vq)) — > as t — > 00. 
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Proof. We refer to Theorem 4.3.3 of |H| for details. Note that u(vq) is the intersection of a 
decreasing collection of nonempty compact sets. Note, in addition, that u(vq) is positive invariant 
by definition and is invariant by compactness argument. 



□ 



Theorem 5.13 Let V be a Liapunov function for t > and, for 

£ := IveV: V{v) = ol , (5.29) 



let K, be the maximal invariant set in S. If the orbit 7(t>o) ^es in a compact set in V, then 
S(t)v — > JC as t — > oo. 

Proof. By definition, V (S(t)va) is a nonincreasing function of t and bounded from below, by 
hypothesis. So, lim V (S(t)v ) = v exists. If y G uj(v ), then V(y) = v and, in view of the fact 

t— >oo 

that S(t)y = y, we have V (S(t)y) = v which implies V(t) = and w(vo) G /C. 

□ 

Now, we apply the invariance principle to the problem at our hand. As we will see, if Bq is a 
sufficient large ball around <po = in the cone C (with the induced topology of H^ p ), the invariant 
set K-k = {^(^o)) u <E B } C £ consists of the union of unstable manifolds for the equilibrium points 

4>o, 4>i, ■ ■ ■ , 4>k, with 4>j(x) = / ipj~{y) dy, provided a is such that 2/ (k + l) 2 < a < 2/k 2 holds for 

Jo 



some k G N. Note that the hypotheses of Theorems [5. 12| and |5.13| hold since the orbits of S(t)vo 



are bounded in by Theorem |5.10| and remain in a compact set of p in view of Theorem [3.2| . 
For this the Sobolev embedding theorem is evoked: W 2,2 (— 7r, 7r) C C l+a (— 7r, 7r) with continuous 
inclusion, so v has a continuous representative in Cl + a which belongs to C 2 +° by Schauder estimates 
(see e.g and references therein). Therefore, any solutions m(£, x) = fjf v(t, y) dy of (p.3|) in C has 
a continuously three-times differentiable representative. 

We thus have 

Theorem 5.14 If a > 2, <p = is globally asymptotically stable solution of ( \3.3j ) in 

C= {u6 Hl p n # e 2 p : m(0) = 0, u > and au xx < 1} . 

If a < 2, the origin is unstable in C and there exits an open dense set U C C of initial conditions 
such that lim u(t; Uq) — > <pi for all Uq GW. 

t— >oo 



Proof. It follows from Theorem |5.13| v (t; •) — > u (v(0; •)) C |^ : V = j> in V as t — > oo. 
But, from «, V (i/j) = iff 



a^" - 2a + 2ip = , (5.30) 
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whose solutions are if) = ip and if)f, j = 1, . . . , k, studied in Section f|. We note that (f>j(x 

/ ip+(y) dy>0 ((f)j(x) < 0) for all x E [-tt, tt] and j > 1, since (V>/)' (0) > ((ipj)' (0) < 0) 
Jo 

Multiplying (|5.30|) by if) and integrating over (— tt, tt), gives 



(c#" + 2^cfe = - |H|, /2 < 
if a > 2. The nonlinear term vanishes since, by integration by parts, 

/*7T 

V>V> 2 dx = -2 / ^V 2 . 



This implies r/> = and proves that S(t)vo — ► as t —>■ oo in V. We quote Theorem 4.3.5 in 
[[H|l for the instability assertion. 

Since the spectrum a(L) of the linearized operator around the equilibrium points (see Theorem 
|5.3j ) lies on the real line, all equilibrium points are hyperbolic, S given in ( |5.29|) is a discrete and 
finite set and 

v = y w s (^) 

holds with W S (V') = {""o G V : S(t)vo — ► ^ as t — > oo}. It is proven in that each stable manifold 
W s {ip) is a C 2 embedded submanifold of V (W s (0) is C 3 submanifold of C) and, if if) is locally 
unstable, than W s (if)) has codimension larger than or equal to 1. Therefore, V, and consequently 
C , can be written as a finite union of open connected sets together with a closed nowhere-dense 
remainder. 

□ 

Finally, we show that, for an open set Vo C V given as before, the maximal invariant set 

K. k = |J W n (V) (5.31) 

where W u (ip) = {v E V : S(t)v Q — > if) as t — ► — oo} is the unstable manifold of if). By Theorem 
5.13| , the orbit v(t; Uq) = S(t)vo exists and remains, by invariance, in for all t E R. Therefore, 



lim v(t,uo) = if) exists and if> E JCk, so Kk C U^f£ Since the converse is also true, the 

t^>CG ^ 

equality (|5.31| ) thus holds. 
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